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Abstract 

The aim of the paper is to answer a long-standing open problem on the 
relationship between NP and BQP. The paper shows that BQP contains 
NP by proposing a BQP quantum algorithm for the MAX-E3-SAT prob¬ 
lem which is a fundamental NP-hard problem. Given an E3-CNF Boolean 
formula, the aim of the MAX-E3-SAT problem is to find the variable as¬ 
signment that maximizes the number of satished clauses. The proposed 
algorithm runs in O(m^) for an E3-CNF Boolean formula with m clauses 
and in the worst case runs in 0{n^) for an E3-CNF Boolean formula with 
n inputs. The proposed algorithm maximizes the set of satisfied clauses 
using a novel iterative partial negation and partial measurement tech¬ 
nique. The algorithm is shown to achieve an arbitrary high probability of 
success of 1 — e for small e > 0 using a polynomial resources. In addition 
to solving the MAX-E3-SAT problem, the proposed algorithm can also be 
used to decide if an E3-CNF Boolean formula is satisfiable or not, which 
is an NP-complete problem, based on the maximum number of satisfied 
clauses. 

Keywords: Quantum Algorithm, MAX-E3-SAT,E3-SAT, Amplitude Am- 
plihcation, BQP, NP-hard, NP-complete. 


1 Introduction 

A long-standing open problem in quantum computing is the relationship be¬ 
tween the classes NP and BQP n a. Decision problems are in NP if yes- 
instances have witnesses that can be checked in polynomial time m- The class 
BQP is the quantum computing analogue of the classical class BPP (bounded 
error probabilistic polynomial) [6]. A problem is in BPP if there is a proba¬ 
bilistic classic algorithm (Turing machine with access to random bits) which 
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makes errors (for either yes or no instances) with probability of given a wrong 
answer at most 1/3. The value 1/3 is arbitrary - all that is required is that the 
value is bounded away from 1/2. By repeated runs, the probability of failure 
can be made exponentially small. The problem class BQP replaces the classical 
algorithm with a quantum algorithm [T] . Thus a decision problem is in BQP if 
there is a quantum algorithm for it with probability of being wrong less than 
1/3. 

The common belief concerning the relationship between NP and BQP was 
that NP is not contained in BQP (e.g see chapter 15 of [8]). However, a recent 
paper by one of the current authors has shown that an NP-hard problem (Graph 
Bisection) can be efficiently solved, with low failure probability, by a quantum 
algorithm m- This implies that NP is in fact contained within BQP, as any 
NP problem can be polynomially reduced to an NP-hard problem. The Graph 
Bisection problem is, perhaps, somewhat obscure and much of the presentation 
of that result involves ensuring the balance of the partition, which detracts 
from the main features of the approach. Gonsequently in the current paper, we 
will directly address the classic Boolean Satisfiability problem (SAT) to show 
precisely how constraints, expressed as Boolean formula, can be encoded into 
quantum algorithm. The constraints are entangled with the superposition of 
all possible truth-value assignments and a probability amplification technique 
applied to amplify the assignment which maximizes the number of satisfied 
clauses. 

In particular, we will focus on MAX-E3-SAT [TT], in which each clause con¬ 
tains exactly three literals, and we will show that our quantum algorithm will 
solve this maximization problem with high probability of success. In particu¬ 
lar, it can then be used to solve the decision problem (with high probability). 
Iterating the process allows the probability of failure to be made exponentially 
small. 

A key fact about the MAX-E3-SAT problem is that random truth assign¬ 
ments will satisfy, in expectation, 7/8 of the clauses. A consequence of the 
PGP Theorem is that this cannot be improved upon (more precisely, there is no 
(7/8 -I- e) approximation algorithm for constant e > 0) unless P=NP |7]. That 
MAX-E3-SAT can be solved (with high probability) in polynomial time by a 
quantum computer is therefore all the more remarkable. 

The result shown in this paper doesn’t contradict with that shown [2 where 
it was shown that BQP does not contain NP relative to a random unitary 
oracle with probability one. This argument does not imply that BQP does not 
contain NP in a non-relativized world which is the novel feature in the proposed 
algorithm, where partial measurement is used in the amplitude amplification 
process instead of the usual unitary amplitude amplification techniques that 
use iterative calls to an oracle to amplify the required solution. 

The aim of the paper is to propose a quantum algorithm for the MAX-E3- 
SAT problem. The algorithm prepares a superposition of all possible variable 
assignments, then the algorithm evaluates the set of clauses using all the pos¬ 
sible variable assignments simultaneously and then amplifies the amplitudes of 
the state(s) that achieve(s) the maximum satisfaction to the set of clauses using 
a novel amplitude amplification technique that applies an iterative partial nega¬ 
tion and partial measurement. The proposed algorithm runs in O(to^) for an 
E3-GNF Boolean formula with m clauses and in the worst case runs in 0{n^) for 
an E3-GNF Boolean formula with n Boolean variables to achieve an arbitrary 
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high probability of success of 1 — e for small e > 0 using a polynomial resources. 

The paper is organized as follows; Section 2 shows the data structures and the 
quantum circuit for encoding an E3-CNF Boolean formula. Section 3 presents 
the proposed algorithm with analysis on time and space requirements. Section 
4 concludes the paper. 

2 Data Structures and Clause Encoding 

An n inputs fc-CNF Boolean formula, 

f{Xo,Xi, . . .,Xn-l) = Co A Cl A ... A Cm-l, (1) 

is a conjunction (AND) of m clauses, each clause cj represented by a disjunction 
(OR) of exactly k < n literals, Cj = [ijg V V ... V such that a literal 

Ij^a in clause Cj with 0 < a < /c — 1 and 0 < j < m — 1 equals to an input variable 
in its true form Xi or its complemented form -^Xi. That is, Ij^a = Xi, where Xi can 
be replaced by Xi or ^Xi such that -icci is the negation of Xi with 0 < i < n — 1. 
The first aim is to decide whether / is satisfiable or not (deciding /). The 
second aim is to find a variable assignment for Xo,Xi,... and Xn-i that satisfies 
/ if it is satisfiable (solve /) or to find a variable assignment that satisfies the 
maximum possible number of clauses if / is unsatisfiable (maximize /). 

The problem of deciding whether a Zc-CNF Boolean formula is satisfiable or 
not is NP-complete and is known as fc-SAT or Efc-SAT problem. The optimiza¬ 
tion problem associated with the Ze-SAT problem to find a variable assignment 
to satisfy a satisfiable fc-CNF formula is NP-hard. If the Zc-CNF is unsatisfiable, 
then the problem of finding a variable assignment to maximize the number of 
satisfied clauses is known as MAX-Ek-SAT problem which is an NP-hard prob¬ 
lem m- The maximum number of clauses for a fc-CNF Boolean formula is 

= 0{n^). Without loss of generality, this paper targets the E3-SAT 

and MAX-E3-SAT problems where fc = 3, so the maximum number of clauses 
m for the MAX-3E-SAT is |n(n - l)(n - 2). 

2.1 Encoding of a Solution 

A candidate solution S to the MAX-E3-SAT problem is a vector of variable 
assignment A = (xq, , a;„_i) € {0,1}", and each vector A is associated 

with a vector of the truth values C{A) = (cq, Ci,..., c^-i) S (0,1}™ of the m 
clauses sorted in order, i.e. S = (A, (7(A)) G (0,1}”+™. The optimal solution 
‘S’max = (Amax, C'(Amax)) IS the solution that contains a vector of variable 
assignment Amax that maximizes the number of I’s in the vector of the truth 
values C{Amax) of the m clauses. For short, the number of I’s in the vector 
of clauses (7, i.e. the number of satisfied clauses, will be referred to as the 1- 
density of C so that the 1-density for a satisfiable formula must be equal to m. 
For example, consider the E3-CNF Formula with n = 3 and to = 4, 



f{xQ,Xi,X2) = Co A Cl A C2 A C3, 


( 2 ) 


3 


where, 


Co = (^cco V -icci V -1X2), 

Cl = (^Xo V XiV X2), .ON 

C2 = (a:o V -^Xi V X2), ^ 

C3 = (xo V XiV X2), 

then a solution to this formula will be encoded as S' = {A,C{A)), where 
A = {xo,xi,X2) and C{A) = (cq,C i,C2,C3). This formula is satisfiable when 
{xo,Xi,X 2 ) = (0, 0,1), (0,1,1), (1,0,1), or (1,1,0), and an instance of an op¬ 
timal solution will be Smax = ((0, 0,1), (1,1,1,1)) with A^ax = (0,0,1) and 
C{Amax) = (1,1,1, !)■ For n > 3 and m > 7, the E3-CNF formula might not be 
satisfied |S] where the 1-density of the C{Amax) vector will give the maximum 
number of satisfied clauses and the order of I’s will show the satisfied clauses 
using the variable assignment Amax- 

For n > 3 and m = |n(n — l)(n — 2), we have the problem instance comp- 
srising all possible clauses. The 1-density of C{A) in this case will be 7m/8 
which is the worst possible case for the 1-density of C. The 3-CNF formula will 
be unsatisfied for an arbitrary variable assignment ^ M- 

2.2 Encoding of a Clause 

An E3-CNF formula with n inputs and m clauses will be encoded as an n -I- m 
inputs/outputs quantum circuit. Every E3-CNF clause c = {Iq \/ li V I 2 ) will be 
encoded using a 4 x 4 quantum gate. The GT^ (4x4 Generalized Toffli) gate 
[S] is the main primitive gate that will be used to encode a clause. The GT'^ 
gate is defined as follows: 

Definition 2.1 (GT^ gate) 

GT^ gate is a reversible gate denoted as, 

(2/O) 2 / 1 : 2/2; font) = GT'^^Xq © ( 5 o: Xi (B Si, X2 (B S2', fin), ( 4 ) 

where Xa,Sa,fin and font G {0:1} with a G {0,1,2}. The GT'^ gate has 4 
inputs: Xq, Xi, X 2 (known as eontrol qubits) and fin (known as target qubit). 
Eaeh eontrol qubit Xa is associated with a condition Sa, such that if Sa = ^ then 
the condition on Xa is satisfied if Xa = 0, i.e. cCq © 1 = -^Xa, and if Sa = 0 then 
the condition on Xa is satisfied if Xa = 1. The GT^ gate has 4 outputs: yo, yi, 
j /2 and font- The operation of the GT^ gate is defined as follows, 

= a;a, fora = {0,1,2}, 

font = fin © ((a^O © © (^l) A {X 2 © ^2)) : 

where © is the XOR logic operation, i.e. the target qubit /„ will be flipped 
if and only if each control qubits Xa satisfies its associated condition Sa- For 
example, font = ^fin for the gate GT'^^Xq © l,a:i,X 2 © 1; fin) if and only if 
xq = 0, xi = 1 and X 2 = 0. 

A GT^ gate with its target qubit, fin, initialized to state |1) can be used to 
encode a clause c = (/q V 2i V I 2 ) using the Boolean algebraic identity, 

c = (2o V 2i V I 2 ) = ((^0 © 1) A {h © 1) A {I 2 © 1)) © 1, (6) 


4 


Xo^ —2 -' 

►-c 

?— ■^0) 




) 

t J, d 


, ; 


1 < f 

, 


* ' ^3) 

rrs\ 

ll \ 



Co; 

h \ 


[J 

_ 

ci; 

GTi GT^ G. 

h 


Figure 1: A quantum circuit using GT'^ gates for the E3-CNF formula shown in 
Eqn. and where • on a control qubit means that the associated condition 
is 1 while o on a control qubit means that the associated condition is 0. 


so that fout = c, where la = x, and x can be replaced by x or ^x such that 
-la; = a: © 1 is the negation of x. That is, 

(a;o,a;i,a: 2 ; c) = GT‘^{lo © 1, © 1, © 1; 1). (7) 

For example, consider the following E3-CNF Boolean formula with n = 4 
and m = 3, 


f{xo,Xi,X2,X3)=CoACiAC2, ( 8 ) 

with 

Co = (a;o V -^Xi V -ia;2), Ci = (-'a:o V Xi V -^x ^), C2 = (a;o V a;2 V ^a;3). (9) 

Apply the Boolean algebraic identity shown in equation ^ on each clause, 
then. 


Co = (xo V ~^Xi V -^X2) = (“'a^o A a;i A 3 : 2 ) = (-'a;o A a;i A a;2) © 1, 

Cl = (- 1 X 0 V Xi V - 1 X 3 ) = -I (xq a ^Xi A X 3 ) = (xo A -'Xi A X 3 ) © 1, (10) 

C 2 = (Xo V X 2 V - 1 X 3 ) = ^ (- 1 X 0 A ^X 2 A X 3 ) = (-'Xo A -'X 2 A X 3 ) © 1, 

then each clause can be encoded using a GT'^ gate as follows, 

(xo,xi,X 2 ;co) = GT^ (xq © l,xi,X 2 , 1 ), 

(xo,xi,X3;ci) = GTf (xo,xi © 1,X3,1), (11) 

(xo,X2,X3;C2) = GT^ (xq © 1,X2 © 1,X3, 1) . 

To construct a quantum circuit for this E3-CNF formula, prepare a quantum 
register with 4 qubits to be loaded with the values of xo, xi, X 2 and X 3 , and add 
3 extra qubits initialized with the quantum state |1) so that GTq uses the first 
extra qubit as the target qubit, GT^ uses the second extra qubit as the target 
qubit, and so on, as shown in figure Let f7 be a quantum circuit on 7 qubits 
defined as U = GT^GTfGT^, then. 


(xo,Xi,X2,X3;Co,Ci,C2) = t7(xo,Xi,X2,X3; 1,1,1). (12) 
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Figure 2: A quantum circuit for the proposed algorithm. 


3 The Algorithm 

Given an E3-CNF formula / with n inputs and m clauses. The proposed algo¬ 
rithm is divided into three stages, the first stage prepares a superposition of all 
possible variable assignments for the n variables. The second stage evaluates 
the m clauses for every variable assignment and stores the truth values of the 
clauses in truth vectors entangled with the corresponding variable assignments 
in the superposition. The third stage amplifies the truth vector of clauses with 
maximum number of satisfied clauses using a partial negation and iterative mea¬ 
surement technique. The proposed algorithm uses (n -I- m -I- 1) qubits during 
the three stages. Each of the first n qubits is initialized to state |0), each of the 
m qubits is initialized to state | 1 ), and an extra auxiliary qubit, denoted |ax) 
is initialized to state |0). The qubit \ax) is an auxiliary qubit that will be used 
during the amplitude amplification technique. The amplitude of the state of 
I ax) entangled with every truth vector of clauses will act as an aggregator for 
the amount of partial negation to be applied on that state based on the 1 -density 
of the entangled truth vector of clauses. The probability of finding |ax) = |1) 
when measured will depend of the accumulation of all partial negations applied 
on the states of that qubit. 

The system is initially as follows, 

|V'o) = |0)"G|l)'"G|0). (13) 


1- Variable Assignments Preparation. To prepare a superposition of all vari¬ 
able assignments of n qubits, apply iJ®" 0 on the n -I- m -I- 1 qubits 

= (FT®" G /®™+i) IV^o) 

^ fc =0 

where H is the Hadamard gate, I is the identity matrix of size 2x2, N = 
2 ”, and = (xq, x*,..., x^_;^) € { 0 , 1 }” is the bit-wise representation 
of an integer fc, for 0 < fc < — 1, that represents a variable assignment 

out of the N possible variable assignments. 
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Figure 3: Quantum circuits for the operator followed by a partial measure¬ 
ment then a negation to reset the auxiliary qubit |ax). 


2- Preparation of the Truth Vectors of Clauses. For every E3-CNF clause 
Cj = {lo V V I 2 ), apply a GT^ gate taking qubit j in the m qubits register 


as the target qubit as shown in Section 2.2 The collection of all GT'^ 


gates applied to evaluate the m clauses is denoted C/ in figure]^ then the 
system is transformed to, 


\^j:^) = {Gf®I)\^P^) 

N-l /I c'l 

k—0 

where Gk = (cq, Ci,..., c^_i) £ {0,1}"* is the truth vector for the m 
clauses associated with variable assignment Ak- 


3- Maximization of the Number of Satisfied Clauses. The aim of this stage 
is to find the state \Gk) that contains the maximum number of |l)s. Such 
a state will be denoted \Gmax)- A modified version of the amplitude 
amplification algorithm shown in m will be used for this purpose. Every 
\Gk) is entangled with the corresponding variable assignment |Afc). Since 
the variable assignment \Ak) will not be involved directly in this stage 
and the corresponding \Gk) will not be modified by any operation, then 
for simplicity the system can be re-written as, 

k=0 

Let dk = (Gk) be the 1-density of the state \Ck)- A solution to the MAX- 
E3-SAT problem is to find the state \Gmax) with dmax = max{dk, 0 < 
k< N -1}. 

The aim is to find \Gmax) when |'03) is measured. To find \Gmax)-, the 
algorithm applies partial negation on the state of |ax) entangled with \Ck) 
based on the 1-density of \Ck), i.e. more I’s in \Gk) gives more negation 
to the state of |ax) entangled with \Ck)- If the number of I’s in \Ck) is 
m, then the entangled state of |ax) will be fully negated. 
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Let X be the Pauli-X gate which is the quantum equivalent to the NOT 
gate. It can be seen as a rotation of the Bloch Sphere around the X-axis 
by TT radians as follows, 


X = 


0 1 
1 0 


(17) 


The partial negation operator V is the root of the X gate and 
can be calculated using diagonalization as follows, 




1+t 1-t 

l-t 1+t 


(18) 


where t = ^—1, and applying V for d times on a qubit is equivalent to 
the operator. 


yd _ 


1 

2 


l + t‘^ 1 - ■ 

1 - 1+t^ ’ 


(19) 


such that if d = m, then P™ = X. To amplify the amplitude of the state 
\Cmax)-, apply the operator on as will be shown later, where 
is an operator on m -I- 1 qubits register that applies V conditionally for m 
times on |ax) based on 1-density of |coCi... Cm-i) as follows (as shown in 
figure [^, 


Mx = Cont_V{cy ax)Cont_Vax)... Cont_V[cm-i] ax), (20) 

where the Cont_V(cj;ax) gate is a 2-qubit controlled gate with control 
qubit \cj) and target qubit \ax). The Cont_V{cj\ax) gate applies V 
conditionally on \ax) if \cj) = |1), so, if d is the 1-density of |coCi... Cm-i) 
then. 


(|coCi...Cm_i) (g) |0)) 


|coCi...c„_i) (g) 


l+t‘^ 

2 


| 0 ) + 


2 



( 21 ) 

and the probabilities of finding the auxiliary qubit \ax) in state |0) or |1) 
when measured is respectively as follows. 


Pr{\ax) = |0)) 
Pr{\ax) = |1)) 


1 +r 


2 

1-t^ 


= cosM^), 


( 22 ) 


To find the state \Cmax) in IV’a)) Ih® proposed algorithm is as shown 
in Algorithm 1 and as shown in figure For simplicity and without 
loss of generality, assume that a single \Cmax) exists in {ips), although 
different variable assignments might be associated with truth vectors with 
maximum 1-density with different order of I’s, but such information is not 
known in advance. 














Algorithm 1 Amplify \Cmax) in IV's) 

1; Let I'i/'r) = 1 ^ 3 ) 

2: for counter = 1 — )■ r do 
3: Apply the operator on \ipr)- 

4: Measure \ax) 

5: if |aa;) = |1) then 

6: Let \ipr) be the system post-measurement of |ax) 

7; Apply X gate on \ax) {to reset jax) to |0) for the next iteration} 

8: else 

9: Let \ipr) = 1 ^ 3 ) and restart the for-loop 

10: end if 

11; end for 

12: Measure the first m qubits in \'tjjr) to read \Craax)- 

13: if \Cmax) = |1)®™ then 

14: The E3-CNF formula is satisfiable 

15: else 

16; The E3-CNF formula is not satisfiable where number of |l)’s in \Cmax) 
represents the maximum number of satisfied clauses in order 

17; end if 

18: Measure the first n qubits in |'02) to read the corresponding variable assign¬ 
ment \A^ax) 


We require that Algorithm 1 finds |ax) = |1) for r times in a row. The 
probability of finding |ax) = |1) after Line:4 in the I'** iteration of the 
for-loop is given by, 

PrW(|ax) = |1)) = ^ ^ sin^ (23) 

fc =0 ^ ^ 

The probability of finding ji/jr) = \Cmax) after Line:4 in the 1^* iteration, 
i.e. r = 1 is given by, 

Pr^^\\i)r) = \Cmax)) = ^Sin^ (24) 

The probability of finding \ax) = |1) after Line:4 in the iteration, is 
given by, 

Af-l 

PrM(|oi) = |l)) = ,yy?-. (25) 

E (ftf) 


The probability of finding \ipr) = IC'maa;) after Line:4 in the iteration, 
i.e. r > 1 is given by. 
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Figure 4: The probability of success for an E3-CNF formula: (xq V Xi V X 2 ) A 
(-<xo V -<xi V -<X 2 ) with n = 3 and m = 2 with dmax = 2, i.e. the 
maximum number of satisfied clauses is 2, where = \Cmax)) = 

0.75, Pr(i)(|ax) = |1)) = 0.875, = IC^a.)) = 0.98, and 

PA^\\ax) = jl)) = 0.99. 


Pr^dV^,) = \Cmax)) 


sin 


2r 


^max ^ 

2m / 


N-1 

Y: Sin2("-1) 



(26) 


To get the highest probability of success for PrdV'r) = \Gmax)), the for- 
loop should be repeated until [Pr^ ’'Hlaa;) = |1)) - PrG){\'il)^) = \Cmax))\ < 
e for small e > 0 as shown in figure This happens when, 


N-l 

H (is) S «. (27) 

k—0,k^max 

and since the Sine function is a decreasing function then for sufficient large 

r, 


k—O^k^max ^ ' 



(28) 


where dnm is the next maximum 1-density less than dmax- The values of 
dmax and dnm are unknown in advance, so let dmax = m be the number of 
satisfied clauses, then in the worst case when dmax = ni) dnm = m—1 and 
m = |n(n — l)(n — 2), the required number of iterations r for e = 10“^ 
and A > 0 can be calculated using the formula. 
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0 < sin 


(m — 1)'K 
2m 


(29) 


then, 



< e, 


\ loff(e) 

21og(sin(i^)) 

= (30) 

21 og(cos( 2 ^)) 

>HW=0{m^), 

where 0 < m < |n(n — l)(n — 2). When m = |n(n — l)(n — 2), then the 
upper bound for the required number of iterations r is O (n®). Assuming 
that a single \Cmax) exists in the superposition will increase the required 
number of iterations, so it is important to notice here that the probability 
of success will not be over-cooked by increasing the required number of 
iteration r similar to the common amplitude amplification techniques. 


3.1 Tuning the Probability of Success 


During the above analysis, two problems might arise during the implementation 
of the proposed algorithm. The first one is to finding \ax) = |1) for r times 
in a row which is a critical issue in the success of the proposed algorithm to 
terminate in polynomial time. The second problem is that the value of dmax is 
not known in advance, where the value of Pr^^\\ax) = |1)) shown in equation 
(231 plays an important role in the success of finding \ax) = |1) in the next 
iterations, this value depends heavily on the 1-density of \Cmax)^ h®- the ratio 


m 

Consider the case of a complete E3-CNF formula where the number of clauses 
is m = |n(n — l)(n — 2) and all the |C'fe)’s are equivalent where anyone can be 
taken as \Cmax)- In this case, each clause Cj will be satisfied by 7 variable 
assignments out of 8 possible variable assignment, then dmax = |nT, for any 
\Ck) CO], so Pr^^\\ax) = |1)) is as follows. 


pA^){\ax) = \l))=sm^<^) 

= sin2(I|) (31) 

= 0.9619. 


This case is a trivial case for the proposed algorithm by setting m = dmax 
in root of X to get a probability of success of certainty after a single 
iteration. Assuming a blind approach where dmax is not known, then this case 
represents the worst case m and iterating the proposed algorithm will not 
amplify the amplitudes after arbitrary number of iterations. For an arbitrary 
F3-CNF formula, the actual probability of success will depend of the 1-density of 
\Cmax), i.e. the ratio . In the following, a tuning of Pr'd-){\ax) = |1)) will be 
shown so that we can find \ax) = |1) after the first iteration with an arbitrary 
higher probability of success close to certainty without a priori knowledge of 

dmax- 
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Figure 5: The probability of success for the E3-CNF formula shown in figure 
where PA^'^{\ax) = |1)) is raised from 0.875 to 0.94 by adding a single temporary 
qubit initialized to state |1), i.e. Umax = 1 where (Itpr} = \Cmax}) = 0.75, 
Pr(i)(|ax) = |1)) = 0.94, = \Cmax)) = 0.98, and PA^){\ax) = |1)) = 

0.99. 


For an arbitrary E3-CNF formula, we could interpret the formula for Pr{\ax) = 
|1)) in equation (231 as the expected value of the function. 


^(x) 



(32) 


where x is the proportion of clauses satisfied by a random truth assignment, 
that is, Pr{\ax) = |1)) = £’[(/)(x)]. The bounds for the probability of finding 
I ax) = |1) in the first iteration is as shown in the following Lemma, 


Lemma 3.1 The probability of finding |ax) = |1) in the first iteration is bounded 
as follows, 


Proof 


0.691 < 1 - — < Pr<-^\\ax) = |1)) < sin( —) < 0.981. 
32 16 


Pr^^^(|ax) = |1)) = P 



dfc7r \ 

2m J 


< E 


sin 


/ dkTT 
\ 2m 


(33) 


(34) 


Since sin(x) is a concave function on 0 < x < 7r/2, it follows from Jensen’s 
inequality that. 


Pr^^H\ax) = \1)) 

V 2m 


= sm 


< 0.981. 


For the lower bound we use, 

sin^((l — x)ttI2) > 1 — 


TT^X^ 


(35) 


(36) 
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which follows from the Taylor series taken around x = 0. Then, 


PrW(|aa:) = |l))>iEf=-o'(l 

q ^2 3 


77^(1 —dk/m) 
4 


^ 1_L ^ _ 
^ ^ 16, 

> 1 _L 37r- 

— ^ 16. 


4Nm‘^ 


y^iV-1 ,2 

T,k=o dkm 


37^ 

Stt^ 

1.6 


= 1 + 
= 1 + 


= 1 - §2 > 0.691. 


4iVm 


4Nm 

iTT^ 

32 


(37) 


To overcome the problem of low probability of finding |aa;) = |1) in the 
first iteration, we can add fj-max temporary qubits initialized to state |1) to the 
register \Ck) as follows, 


I CqCi ...Cttj—1) ^ IcqCi . . . Cjyi — \CjyiCyji+1 • • ■ — ^ "> 

SO that the extended number of clauses rriext will be mext = m + ^max a-nd 
V = will be used instead of 1^ = in the operator, then the 

density of I’s will be s ^ ^ probability of success Pr^ax to find 

I ax) = |1) after the first iteration of the for-loop in Algorithm 1, 


PrW (|«^) = |1)) = Na^sin^ / "(> + Mn.ax) \ ^ 

y 2 (m + ^max) J 

then the required number of temporary qubits Umax is calculated as follows, 

Mmax ^ 'kn 

where w = |sin“^ (VPpnax) and Na'^ = 1. For example, if Pr^ax = 0.99, 
then (|ax) = |1)) will be in the neighborhood of 99% as shown in figure 
To conclude, the problem of low 1-density of \Cmax) can be solved with 
a polynomial increase in the number of qubits to get the solution \Cmax) in 
O (TOg„.j) = O (n®) iterations with arbitrary high probability Pvmax < 1 to 
terminate in poly-time, i.e. to read |ax) = |1) for r times in a row. 



4 Conclusion 

Given an E3-CNF Boolean formula with n inputs, the paper showed that BQP 
contains NP in a non-relativized world by proposing a BQP quantum algorithm 
to solve the MAX-E3-SAT problem with m clauses. The proposed algorithm 
encoded every clause as a GT^ gate where 0(n -I- m) qubits are used. The 
algorithm is divided into three stages; the first stage prepares a superposition of 
all possible variable assignments. In the second stage, the algorithm evaluates 
the set of clauses for all possible variable assignments using a quantum circuit 
composed of GT"^ gates so that each variables assignment is entangled with a 
truth vector of clauses evaluated according to that variables assignment. In 
the third stage, the algorithm amplified the amplitudes of the truth vector 
of clauses that achieves the maximum satisfaction to the set of clauses using 
an amplitude amplification technique that applies an iterative partial negation 
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where partial negation is applied to the state of an auxiliary qubit entangled 
with the truth vector of clauses based on the number of satisfied clauses, i.e. 
more satisfied clauses implies more negation to entangled state of the auxiliary 
qubit. A partial measurement on the auxiliary qubit is then used to amplify the 
set of clauses with more negation. The third stage requires 0{rn?) iterations 
and in the worst case requires 0{n^) iterations. It was shown that the proposed 
algorithm achieves an arbitrary high probability of success of 1 — e for small 
e > 0 using a polynomial increase in the resources by adding dummy clauses 
with predefined values to give more negation to the best truth vector of clauses. 

In the same manner, the proposed algorithm can also be used decide if a 
given E3-CNF Boolean formula is satisfiable or not by checking the truth vector 
of clauses, if the m clauses are satisfied then the E3-CNF Boolean formula is 
satisfiable, if not, then the proposed algorithm gives the maximum number of 
satisfied clauses with the corresponding variable assignment. 

The proposed algorithm can easily be extended in a trivial way to solve/decide 
an Efc-CNF Boolean formula by encoding any E/c-CNF clause as a GT^+^ gate 
where it can be shown that the algorithm will require 0{n^^) iterations. 
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